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ABSTRACT 

Most pulsars observed by the Fermi LAT have 7-ray luminosities scaling with spindown power E as L-^ 
{E ■ lO'^'^erg/s)^/^. However, there exist one detection and several upper limits an order of magnitude or more 
fainter than this trend. We describe these 'sub-luminous' 7-ray pulsars, and discuss the case for this being an 
orientation effect. Of the 12 known young radio pulsars with E > lO'^^ergs"^ and d < 2kpc several are 
substantially sub-luminous. The limited available geometrical constraints favor aligned geometries for these 
pulsars, although no one case for alignment is compelling. In this scenario GeV emission detected from such 
sub-luminous pulsars can be due to a lower altitude, lower-power accelerator gap. 
Subject headings: gamma rays: stars - pulsars: general 



1. EVITRODUCTION 

The Large Area Telescope (LAT) on the Fe rmi satellite 
has now detected over 75 spin-powered pulsars d Abdo et al] 
|2bl0a; Romani 2011). Among the w 50 non-recycled en- 
ergetic pulsars there is a clear trend for 7-ray 'efficiency' to 
increase with decreasing spin-down power E, giving a heuris- 
tic 7-ray luminosity 
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This is a natural result for models where the emission is 
produced by a Goldreich-Julian current of charges passing 
through a characteristic potential drop (Harding 1981; Aroni 
120061) . Of course, energy conservation limits < E, and 
as E decreases, the star is unable to maintain the potential 
drop, leading to a 'death zone' below E « 10'^'' — IG'^^crg s^^ 
where th is process st a rts to t urn off. This is portrayed in fig- 
ure 5 of lAbdo et aLl (l2010al) . where most energetic pulsars 
lie between Eq ( 1 ) and unit efficiency. Only two young pul- 
sars in that plot lie significantly below the L^,heu line: PSR 
J0205-I-6449, where a small inferred distance places it just be- 
low this value, and PSR J0659-I-1414 (to be discussed in this 
paper) which is ^ 20 x less luminous. Thus, independent of 
its physical validity, Eq. (1) forms an effective lower lumi- 
nosity envelope to the bulk of the observed pulsar sample. 

Estimates of L-^ suffer two complications. The first is the 
source distance; for most LAT pulsars we have only distance 
estimates based on the pulsar Dis persion Measure (DM). DM 
modeling jCordes & Lazidl2002l hereafter CL02) is believed 
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to provide statistically useful estimates of pulsar distances, 
with a scatter of « 30% about independent distance estimates, 
although typical errors for nearby pulsars may be as large as 
60% (see Deller 2009). DM distances ar e certainly not reli- 
able for individual objects, and it appears (A bdo et aLll2010ah 
that they may be especially poor for the young, energetic LAT 
pulsars. This is likely since the sample is nearby and asso- 
ciated with regions of active star formation where the excess 
ionized gas may significantly perturb the dispersion measures. 
About a third of the LAT pulsars are found directly in the 7- 
ray data through so-called 'blind' searches dAbdo et al]|2009l ; 
ISaz Parkinson et al.l2010l) ; most of these lack radio detections 
and so do not even have DM distance estimates. The second 
complication is the conversion from the observed energy flux 
Fe along the Earth line-of-sight to the true sky averaged lu- 
minosity 

= AttMEeD^. (2) 

IWatters et al.1 (l2009h and Ro mani & WattersI ((201^, RWIO) 
have estimated 'flux conversion factors' /n for this correc- 
tion for a variety of pulsar models and viewing geometries. 
For most of the observed pulsars, /q should be in the range 
0.7 — 1.3, altho ugh some lower E puls ars, especially 7- 
selected objects (Wa tters & Roman"! 2011), may have as 
small as 0.1 for 'outer gap' (OG) geometries. 

However, there are a handful of pulsars whose observed lu- 
minosity or limit fall an order of magnitude or more below 
L-f^heu- In spite of the uncertainties just discussed we can 
make a case that they are truly sub-luminous. There are three 
possible interpretations. The first is that the 7-ray radiation 
is beamed away from the Earth line-of-sight (or equivalently 
fn > 10). The second is that some particular physical prop- 
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Fig. 1 . — The spin-down-luminosity plane for energetic pulsars, with the heuristic luminosity trend, which saturates somewhere in the 'death zone' (shaded). 
Unpulsed (DC) E > 0.1 GeV luminosities or limits are plotted, assuming /n = 1 (the RWIO /n = 0.13 point for Geminga is also shown). Left: Objects 
with parallax distance measurements. The 95% error bars for LAT-detected objects (circles) include the flux imprecision, but are dominated by the parallax 
uncertainty. For the 95% upper limits (triangles), the en'or flags represent the parallax uncertainty. Right: Objects with DM-estimated distances. Circles: The 
DC luminosities for radio pulsars with LAT pulse detections. Squares: LAT DC detections (this paper). Triangles: DC upper limits. All eiTor bars include an 
assumed 30% DM distance uncertainty. For a few of the fainter LAT detections flux uncertainties contribute significantly. For PSR J0745— 5353 the luminosity 
at the DM-estimated distance is 10 X lower than the point shown. 



erty of the pulsar prevents them from producing the bright 
high-altitude 7-ray emission typical of other energetic pul- 
sars. The third is that the DM distance is especially poor and 
the pulsar is much more distant than estimated. We test here 
the first possibility, that 7-ray beaming explains the low ob- 
served fluxes of some nearby energetic pulsars. We also com- 
ment briefly on the possibility that objects with detected lu- 
minosities ^ L~fjieu may be probing an emission component 
different to the powerful high-altitude gap emission which ap- 
parently dominates the bulk of the LAT-detected pulsars. 

2. THE SUB-LUMINOUS PULSAR CANDIDATES 

To find sub-luminous pulsars, we measure the DC (un- 
pulsed) flux at the positions of nearby (d < 2 kpc), energetic 
(E > 10'^'' erg s^^) non-r ecycled radio pulsars se lected from 
the ATNF pulsar catalog (iManchester et al.ll2005l) . There are 
12 such objects (Table 1 also includes two comparison ob- 
jects). Since the LAT has detected several pulsars, especially 
millisecond pulsars, in the ^ lO'^^^'^'^erg boundary of the 
'death zone' we also consider the well-studied nearby E — 
jQ33.6gj.gg-i pjjjg^ J1932H-1059 (B1929H-10), which 
has a low LAT flux limit. Finally, for comparison we include 
Geminga (J0633-I-1746), a nearby 7-selected pulsar with an 
HST parallax measurement. We should note that this distance 
cut-off is somewhat arbitrary; for example PSR J1747— 2958 
with a CL02 distance of 2.01 kpc is a LAT pulsed detection. 

To measure the unpulsed fluxes, we use 24 months of LAT 
data (Aug 4 2008 - Aug 4 2010) and the P6_V11 instrument 
response function, a refinement to previous analyses reflect- 
ing improved understanding of the point spread function and 
effective area ( Abdo et al. 201 1). 'Diffuse-Class' events were 
selected from good runs with rocking angle < 52°, recon- 
structed energies —0.75 < Log (i?..^ /GeV) < 2, and a recon- 



structed zenith angle < 100°. The list of point sources used 
in the background model is drawn from a preliminary version 
of the two-year Fermi catalog. The analysis used an updated 
version of the model for the diffuse background - Galactic, ex- 
tragalactic, and residual cosmic rays - that is being prepared 
for publication by the LAT team. Like the model used for the 
IFGL catalog (Abdo et al. 2010b) it is based on fitting tem- 
plates for the diffuse emission to the LAT data. 

For each pulsar we assume an exponentially cutoff spec- 
ti-um dN/dE = Na (E/GcY)-^ cxp{-E/Ec). For the 
bright LAT-detected pulsars (marked ^ in the Table) we al- 
low Ec and F to vary in the fits; the results are consistent 
with parameters quoted in Abdo et al. (2010a). For the other 
pulsars we set these parameters to values determined from an 
empirical fit to detected LAT pulsai's (RWIO): F = -4.1 + 
0.156 logio E and Ec/GeV = -0.45 + 0.71 logm Blc, with 
i?LC the magnetic field measured at the pulsar's light cylin- 
der. We evaluate the likelihood for A^o 3l the known puls ar po- 
sition using 'pointlike', a binned likelihood analysis tool dKenj 
IMS, and using a Bayesian approach with a uniform prior we 
integrate the likelihood to 97.5% to obtain a 2a upper limit on 
the flux. For sources with apparent DC emission, we deter- 
mine the corresponding 9 5% range for the me asured A^o- For 
comparison with results in lAbdo et alJ(l2010al) these measure- 
ments and upper limits are then converted to _E > 0.1 GeV 
fluxes using the model spectra. The uncertainties reported for 
the measured fluxes are statistical only; additional systematic 
error arises from uncertainty in the effective area of the LAT 
(about 20% below 1 GeV, 10% at 1-10 GeV, and 30% above 
10 GeV) and the structure of the diffuse background. Sys- 
tematic uncertainties in the upper limits stem primarily from 
uncertainty in the background model and are comparable in 
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magnitude to those associated with the assumed beaming fac- 
tor discussed below. 

Our study gives six candidate sub-luminous pulsars 
(marked in bold). Three in the uniform sample (plus PSR 
J1932H-1059 = B1929H-10) have parallax distance measure- 
ments. These are particularly important as the parallax con- 
straints control a major factor in the luminosity uncertainty, 
allowing us to probe the effects of beaming geometry and gap 
emissivity. For the others we must rely at present on the DM 
distance estimates. These pulsars are displayed in Table 1 
and Figure 1 . Figure 1 also shows several other nearby non- 
recycled LAT-detected pulsars, highlighting the separation of 
our sub-luminous set from this sample. For this figure we 
have assumed — 1 for all pulsars. The plotted luminosity 
errors are dominated by the distance uncertainties, but do in- 
clude the statistical flux errors. Of course, systematic errors 
and non-unity may add additional uncertainty for individ- 
ual pulsars. 

3. EXTERNAL ANGLE CONSTRAINTS 

For simple dipole models (e.g. the OG model) the pulse 
profile and the expected radiation on the Earth line-of-sight 
are determined by the magnetic inclination angle a and the 
viewing angle If these angles are known, we can predict 
7-ray pulse profiles and fluxes for specific models and correct 
observations to the true L^. Unfortunately these are poorly 
known in many cases. 

3.1. Radio Polarization Data 

The sub-luminous candidates treated here are known ra- 
dio pulsars, so the magnetic impact angle /3 = C — a is be- 
lieved to be small. In the context of the rotating vector model 
(iRadhakrishnan and Cookdl 19691) radio polarization data can 
constrain the viewing angles. In most cases, the small range of 
phase illuminated by the radio pulse allows only an estimate 
of the magnetic impact angle 

/3 = C - a « sin-i[sina/(d5'/d,^)„,ax] (3) 

where the maximum rate of the polarization position angle 
(PA) sweep ^'((/') occurs at ^d^j.max, near the closest approach 
to the magnetic axis. Here the sign of the sweep is meaning- 
ful, determining whether the line of sight is closer to or far- 
ther from the positive rotation axis than the observed magnetic 
pole (at inclination a). Occasionally, when the radio pulse is 
very broad or when the pulse profile presents an inter-pulse, 
the radio polarization can make meaningful estimates of both 
a and from fits to the full polarization sweep 

tan(* + ^-o) = — ^ ' , . (4) 

smC^ cosa — cosC sma cos(0 — (po) 

wh ere the polar ization has the absolute position angle $0 at 
4>o- iKeith et aP (2010) have recently presented several exam- 
ples of constraining fits of Eg. (4) to high quality polariza- 
tion data. As described by lEverett & Weisberg ( 2001), while 
nearly all authors fit to Eqs. (3) and (4), given the standard 
astronomical convention of position angle measurement (in- 
creasing N through E) these equations are inconsistent with 
pulsar angles increasing from the positive spin axis (the 'RVM 
convention problem' ). To be consistent, one must actually use 
aEWOi = TT - urvm and Pewoi = -Prvm- Usually this 
correction is only a formality, but as fits to the 7-ray emission 
improve, including details of sweep-back and magnetospheric 
currents, the signs can be important. Thus in the figures and 



discussio n to follow, we convert all 'RVM' -fit angles to the 
consistent lEverett & Weisberg (2001) convention; we encour- 
age future workers to do the same. 

Other phenomenological constraints may be extracted from 
the radio data. For example, radio emission is generally be- 
lieved to be produced within the 'open zone' above the polar 
cap. For a static aligned dipole the half opening angle covered 
by this radio beam is 

P = VKc (5) 

radians for modest emission altitudes h^c = 2TTh/Pc. If the 
observed radio pulse fills this cone we can write Hlc in terms 
of the pulse width W ~ 2p 

= 4/9 Acos^ [cosa cosC + sinasin^cos(VF/2)]; (6) 

if the radio emission does not fill the open zone this provides 
a lower limit for the emission height. It has also been shown 
that, due to a combination of field line sweep back and aberra- 
tion, the phase of the center of the radio pulse should lead 
the phase of the max PA sweep rate by 

A(/.pA « 4/iLc (7) 

(iBlaskiewicz. Cordes & Wassermaniri99U iDyks 2008. eg.). 
Observationally we identify (j>j with the mid-point of the pulse 
at 10% of its peak and 0pa is identified with 0o in an RVM 
fit. The true phase of minimum magnetic angle is between 0/ 
and 0dV).max- These expressions assume simple static dipoles 
and low altitudes. We have checked against detailed numeri- 
cal simulations of swept-back dipole magnetospheres and find 
that the actual pulse intensity center and phase of maximum 
PA sweep are both sensitive to details of the magnetic field 
structure, especially conditions at the light cyl inder that de- 
fine the edge of the open zone (iCraig et alJl20i ll) . These dif- 
ferences are modest at h^c < 0.05. For objects indicating 
higher altitude radio emission detailed comparison with the 
numerical results can be important. 

In practice, radio pulse p rofiles may represent 'patchy ' illu- 
mination of the radio zone dLvne and Ma nchester 1988), even 
for these young pulsars. This complicates our estimates of W 
and (f)j. In Table 1, we list both Wio, the full width of the ra- 
dio pulse at 10% of the peak intensity and Wi, an estimate of 
the pulse width at 1 % of the peak. These measurements were 
made on archival 1.4 GHz profiles (see Table references). The 
Wi estimate is necessarily approximate, especially for the 
lower S/N pulse profiles. At such low flux levels, extended 
pulsed emission may be generated by interstellar scattering 
tails, weak emission components unassociated with the main 
dipole cap or even non-linearities in the measurement system. 
Nevertheless, for at least a few of these pulsars, this broader 
width captures weak components of the pulse coming from 
the principal emission zone. Further, in some cases, the as- 
sumption of pure dipole geometry and even the identification 
of the radio beam with the open zone are suspect. However, 
despite all of these caveats, these radio measurements do pro- 
vide some phenomenological constraints on the range of al- 
lowable a and even when values for individual pulsars are 
suspect. 

3.2. PWN Torus fits 

When the pulsar wind momentum is equatorially concen- 
trated, one may observe a 'torus' of emission at the spin equa- 
tor, as for the Crab and Vela pulsars. The most useful exam- 
ples are found in Chandra (CXO) X-ray images of PWN tori. 
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TABLE 1 

Young Local, Energetic Radio Pulsars: DC fluxes and Geometry Constraints 



Name 




log(rc) 

[y] 


[kpc] 


-fE,>0.1GGV 
[lO^^^ergcm^ s^i] 


[deg] 


Wi 
[deg] 


[degf 




Blc 
[kG] 


Ref. 

e 


J0358+5413 


34.65 


5.75 


1.04+0-?! 

— 0.16 


< 12.8 


38.9 


56 


13.5 


0.059 


2.1 


1 


J0534+2200'' 


38.66 


2.9&f 


2.0 


1828.^21. 










980. 




J0538+2817 


34.69 


4.60-f 


^•'^'-'-0.16 


< 4.5 


55.4 


92 


35.0 


0.153 


2.3 


2 


10633+1746°-'' 


34.52 


5.53 


25+°-12 
'^•^''-0.06 


4340.+^^' 










1.2 




J0659+1414'' 


34.58 


5.05 


n 9SS+0 033 




31.3 


54 


13.9 


0.061 


0.8 


3 


J0745-5353 


34.04 


6.10 


0.25 


< 0.98 


34.9 


61 


4.0 


0.017 


0.8 


unpub. 


J0834-4159 


34.99 


5.64 


1.66 


< 5.6 


24.0 








3.9 




J0835-4510'' 


36.84 


4.05 




9466.13: 


16.9 


27 


4.2 


0.018 


45. 


4 


J0857-4424 


34.42 


5.35 


1.94 


iQ 1+6.6 


18.7 








0.8 




J1057-5226'' 


34.48 


5.73 


0.72 


294+11 


31.0 








1.3 




J1722-3712 


34.52 


5.53 


1.85 




13.3 


22 


12.0 


0.052 


1.2 


unpub. 


J1740+1000 


35.36 


5.06 


1.24 


< 2.9 


42.5 


71 


28.0 


0.122 


4.7 


unpub. 


J1932+1059° 


33.59 


6.49 


361+0 010 


< 2.9 


19.1 


55 


4.5 


0.020 


0.4 


4 


J2043+2740'' 


34.75 


4.08-^ 


1.80 


1 9 7+3.6 

'-3.0 


16.6 


32 


8.0 


0.035 


3.7 
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" Comparison pulsars, not members of the uniform radio-loud, E > lO^^erg s <2 kpc set. 
Fermi LAT pulsed detection. 

Parallax distances from the ATNF data base. CL02 DM distances (w/o eiTors). Classical Crab kinematic distance. 

Aijipyi measured from the Wio pulse center. 

° Reference for the polarization profile used to fit the pulse widths, offsets and polarization sweeps: 1= Go uld & Lvnel < 19981) . 
2=von Hoensbroech & Xilouris 1 1997), 3=:Everett & Weisbera (2001), 4=Johnston, et al. (2005), 5=Noutsos et al. (20TJ) 
^ Age of the associated supernova remnant; is larger, implying initial sp 



where synchrotron emission is produced in the mildly rela- 
tivistic flow downstream from the termination shock. Doppler 
boosting allows one to distinguish the 'front' and 'back' sides 
and so fi tting can measu re the spin axis inclination to the line 
of sight dNg & Romani|[2008l) . The images are not, however, 
sensitive to the sign of the spin, so a torus fit cannot distin- 
guish between C and C — n — Occasionally, symmetric 
jets also allow C, estimates. These fits provide relatively ro- 
bust model-independent constraints on C,, largely orthogonal 
to the RVM measurements, so that the combination provides 
a good picture of the pulsar geometry. 

Unfortunately, to date no strongly 'sub-luminous' pulsar 
has an X-ray torus measurement, since these are typically 
available only for very young t < 10^'^y pulsars. How- 
ever, PSR J1930H-1852 (unseen by the LAT) has a fit angle 
C, = 33 ± 3° su ggesting that its OG emission should not be 
visible dNg & Ro mani 2008); as the LAT exposure increases 
this can eventually be a useful comparison. Also there is some 
hope of obtaining C, estimates for older (even millisecond) pul- 
sars from fits to the geometry of Ha bow shocks; in some 
cases (eg. Romani et al. 2010, J1741— 2054) these show clear 
signs of equatorially concentrated momentum flux and thus 
opportunities to constrain C,. 

4. INDIVIDUAL OBJECTS 

We can combine these various geometry constraints to re- 
strict the viable location of a pulsar in the (a, Q plane and 
for this location compare with the predictions of the various 
7-ray emission models. As an example, we show in Figure 
2 the constraints for the bright, well-studied Vela pulsar PSR 
J0835-4510 (B083 3-45). The RVM model was or iginally 
developed for Vela ( Radhakrishnan and Cookdll969h so it is 
not surprising that the high S/N Vela data provides good con- 
straints on 13. As usual a is not well determined. However, 
there is a go od measurement of C from CXO images of the 
X-ray torus dNg & Romanil 120081) . Finally the pulse width 
constraints are shown. The background gray-scale shows the 
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Fig. 2. — The spin geometry plane for Vela, with various observational 
constraints. The background gray scale shows the goodness of fit of the LAT 
light curve to a basic OG model, with dark colors better agreement. 



'goodness of fit' of the LAT Vela profile to an Outer Gap light 
curve computed for the particular a and C,, assuming a current- 
free retarded dipole field structure (RWIO). 

The contours marked RVM are quite crowded, since the 
high S/N Vela data provide very strong (3 constraints. These 
cross the C, constraints from the X-ray torus and the two 
constraints select solutions at (a, Q = (56°, 63.5°) and 
(109°, 116.5°). For the estimated A^p^ and the result- 
ing h^c ~ 0.018, the pulse width Wio is easily consistent 
with these solutions. They are not consistent with the Wi 
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Fig. 3. — The spin geometry plane for sub-luminous pulsars with parallax distances. The backgrounds show the generic locations providing sharp OG pulses, 
except for PSR J0659+1414, where the background shows the allowed fits to the observed LAT pulses, including lower altitude (eg. 'Two Pole Caustic' TPC) 
emission. Three green contours show the loci of best RVM matches, while the bold and narrow magenta curves showed the regions allowed by emission from the 
static dipole open zone for our estimated Hlc- For PSRs J0659+1414 (B0656+14) and J1932+1059 (B1929+10) the RVM fits of Everett & Weisberg (2001) are 
indicated. For PSR J0538+2817 the fits imply large hi^c requiring a numerical magnetosphere model. The fits to the polarization geometry using such models 
are shown by the dashed (red) contours. 
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width, but such a wider pulse could be easily accommodated 
if this faint emission comes from slightly higher altitudes. 
Indeed, evidence of variable pulse co mponents and micro- 
structure in the Vela pulse wings (J ohnston, et alJi200 1) sug- 
gests such multi-altitude emission. The darker gray scales in 
the background show regions with good fits to the LAT light 
curve. Evidently, the best regions in this model are not at the 
RVMh-PWN preferred positions. As shown in RWIO, model 
perturbations such as magnetospheric currents can shift the lo- 
cations of the best fits; for example OG models with currents 
can produce reasonable agreement with the (109°, 116.5°) 
solution. However other models remain viable. For exam- 
ple the 'Separatrix Layer' (SL) model of iBai & Spitkovskvl 
d201Q,) produces a good match at the (56°, 63.5°) solution. 

4.1. Pulsars with Parallaxes 

In Figure |3] we show the constraints for the sub-luminous 
candidates having parallax measurements. The background 
gray scales show the region where a narrow OG, with gap 
widths w « L^,heu/E, produces sharp caustic pulses. The 
gray levels indicate goodness of fit for a generic single 7-ray 
pulse. Of course with an actual LAT detection the detailed 
light curve and phase produce much more detailed constraints 
within this envelope (see Figure|2]i. Note that sharp OG pulses 
are not expected (white background) near the spin poles un- 
less the pulsar is an orthogonal rotator. In contrast, lower al- 
titude (e.g. TPC) models produce emission at small /3 all the 
way to the poles and good pulse matches are expected in the 
OG blank zones, if such low altitude gaps are active. 

PSR J0358H-5413 (B0355H-54) is a bright radio pulsar for 
which RVM fits provide good (3 constraints and a preference 
for a > 110°. At the estimated A0p^, the Wiq pulse width 
does not significantly tighten this bound, but if we include 
the wider Wi width the constraints tighten with a > 130° 
and C > 140°. This is also the region of best RVM fits. 
Thus while the best-fit geometry localizes to {a, C,) where 
OG emission would not be visible, acceptable solutions in- 
clude the range that could be consistent with visible high alti- 
tude 7-rays. Since our present flux bound only restricts us to 
< Ljjieu/i we regard this as a plausible sub-luminous pulsar, 
but not a strong case. Increased LAT exposure and improved 
geometrical constraints are needed to make this definitive. 

PSR J0538H-2817 in contrast has only rather poor RVM 
constraints. However, the large pulse width does help re- 
strict the range of viable solutions, even though /ilc is high. 
Here the best combined radio data seem to prefer small a, 
small In fact if the Wi constraint is included we conclude 
a < 35°, C < 50° so we would not expect to see OG emis- 
sion. However, using only the Wio width a wide range of 
a is allowed towards the edges of the RVM and pulse width 
contou rs. CXO in fact show s a small X-ray PWN around the 
pulsar (iRomani & Ngl 120031: iNg et al.ll2007h . Unfortunately, 
the emission is too faint and compact to provide a good ( 
measurement, although the existence of opposing jet- like fea- 
tures suggest ss 90°. Interestingly, iKramer. et al.l (l2003h 
find an RVM fit giving a ss 85° and /3 w -2°, which would 
be consistent with jets viewed near side-on; however this so- 
lution is well outside the RVM-allowed region in Figure [3] It 
is worth noting that with t his pulsar's l arge h r.c = 0.14 the 
detailed field treatment of ICraig e t al.l (1201 ll) can be impor- 
tant. The dashed (red) contours show the constraints from fits 
to these numerical models. The best fit altitude depends on a 
and C, ranging from /ilc — 0.13 — 0.17. The contours show 



best solutions to these numerical models, marginalized over 
fiLc- Good solutions are found for a < 40°, acceptable solu- 
tions are at a < 80°. Again we must conclude that alignment 
is preferred, but nearly orthogonal rotators are not excluded. 
This is a case where improved radio observations with higher 
S/N can substantially improve the polarization modeling con- 
straints. 

The next object, PSR J0659H-1414 (B0656H-14) is the 
archetype sub-luminous pulsar. This object was included 
in the very careful polarization study of lEverett & WeisbergI 
(l200lh . who fi nd a = 29 ± 23° and = 8.9 ±6.1° (blue el- 
lipse). Earher. H^ne and Manchested (119881) found a k. 8.2° 
and (3 k, 8.2°, so there has been some consensus that this is 
an aligned pulsar. However, we and Weltevrede et al. (2Q\(% 
find less strong constraints on a even though we are fitting 
the same 1 .4 GHz data from long Arecibo integrations used in 
lEverett & WeisbergI (|200T1) . In our case, the best fits slightly 
prefer small a < 80°, although all values are allowed. How- 
ever, PSR J0659H-1414 has weak but well measured emission 
extending well beyond the peak of the pulse, giving a large 
Wi ~ 54°. Including this constraint does indeed prefer near 
alignment. Our best fits are in fact for a < 35° although the 
large a > 130° fits are nearly as good. As the background 
to the constraints in this panel we show the gray scale good- 
ness of fit for the LAT data compared with a TPC model. This 
model includes emission from below the null charge surface 
and so predicts gamma-ray detections all the way to the spin 
axis. There are in fact two broad zones of reasonable TPC 
model fits consistent with the radio constraints near a w 40° 
and a « 140°. In contrast the OG model for this pulsar's E 
has no 7-ray emission for solutions at a < 50° or a > 130°. 
Thus, if we adopt both the RVM and pulse width constraints, 
a large fraction of the remaining phase space is incompatible 
with 7-ray emission from an OG model and a lower altitude 
(eg. TPC-type) component is preferred. If we adopt the WI 
constraint or the lEverett & WeisbergI (120011) a value this be- 
comes a strong conclusion. However with the less strict angle 
constraints inferred in this paper, OG exclusion is suggested, 
not required. 

The final parallax candidate is PSR J1932H-1059 
(B1929H-10), which has lower E than the cuts for our 
uniform sample. This pulsar has been the subject of many 
polariz ation studies, summarized in Everett & Weisberg 
(1200 lb . who find a = 35.97 ± 0.95°, ( = 61.52 ± 1.3° from 
a fit excluding the main pulse, with the high S/N leading to 
very small statistical errors (blue circle). For consistency 
with the other objects in this paper, we have fit the main 
peak width and position angle sweep. Our RVM fit prefers 
a < 60°, close to the EWOl value, while allowing all a, as 
usual. However, given the large pulse width and low /ilc 
of Table 1, nearly aligned rotators are preferred. For the 
Wio width we infer a < 20°, while the Wi width implies 
a < 15°. The situation for this pulsar is not clear; these 
pulse width constraints are not consistent with the EWOl 
fit. Moreover, this pulsar has a widely separated faint pulse 
component which would be identified as an interpulse for 
orthogonal solutions, but for more aligned solutions suggests 
a very wide pulse profile. Such emission can only come 
from the open zone for very large h^c, which we do not 
model here. Nevertheless, for either the EWOl solution or the 
RVM/pulse width constraints, high altitude 7-ray emission is 
not expected to be visible at Earth for this small E pulsar At 
first sight this would seem to provide a strong confirmation 
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Fig. 4. — The spin geometi'y plane for two pulsai's without parallax distances, showing the radio polarization and pulse width constraints. For geometries away 
from the gray background, the sources are not expected to have strong outer magnetosphere 7-ray pulses. The PSR J1740+1000 data suggest large /i^c requiring 
numerical modeling; the locus of best fits for these models is shown by the dashed (red) contours. 



of the connection between alignment and low 7-ray flux. 
However, with such a low E, this pulsar is in the 'death zone' 
where powerful 7-ray gap emission may have turned off. 



4.2. Sub-Luminous Candidates without Parallaxes 

Since the objects with precise distances do not yet provide a 
definitive test of the nature of sub-luminous pulsars, we check 
other pulsars for which the LAT provides relatively low lu- 
minosities at their estimated distances. For example, PSR 
J0745-5353 (B0743-53) is assigned a distance of 0.25kpc 
in the CL02 model, because it is superimposed on HII emis- 
sion associated with the Gum nebula. At this distance the flux 
is > 450 X less than that expected f rom L^j,.pji- How e ver, at 
the 7.1 kpc distance implied by the iTavlor & CordesI (119931) 
model the upper limit is not constraining (the pulsar does 
remain sub-luminous for distances as large as 2 kpc). The 
combined RVM and pulse width constraints imply a > 150° 
(H/10) or a > 160° (Wi). Clearly, these constraints indicate 
an anti-aligned rotator such that only 7-ray emission from 
below the null charge surface (r < tnc) should be visible. 
Unfortunately the highly uncertain distance prevents us from 
drawing strong conclusions from the absence of flux from this 
pulsar A parallax distance measurement (or at least a lower 
limit) would be of particular value; if a low distance is con- 
firmed it provides a sharp test of the 7-ray beaming geometry. 

Our final pulsar is PSR J1740H-1000. This pulsar is of inter- 
est since it is young and is located well off the Galactic plane. 
This makes for a very clean LAT flux limit and a relatively 
robust DM distance. We have performed RVM fits using po- 
larization data collected at Nangay with the BON pulsar back- 
end (see Theureau et al. 201 1). The fits allow a large a range. 



but prefer values near 120°. However, such an orthogonal 
rotator is very difficult to accommodate with the very wide 
observed pulse, even for the relatively large hic = 0.122 in- 
ferred here from A(j)pA- The standard pulse width constraints 
are shown in Figure 4; to accommodate the Wia width implies 
a < 70° or > 120°, while the Wi width requires a < 50° 
or > 140°. Interestingly, if we fit the numerical magneto- 
sphere models, the large a solution becomes preferable and 
in fact shifts to slightly higher values. The red (dashed) con- 
tours show the numerical model fit including ^I* points within 
WiQ. For numerical fits placing points out to the Wi pulse 
width in the open zone, the allowed region moves inside the 
magenta Wi contour and we find a w 30° or 150°, where 
one would expect no high altitude 7-ray emission at this E. 
However, if we only consider the data within the Wia pulse 
width an appreciable region including OG-type emission is 
allowed. Thus the geometry of this pulsar is not yet suffi- 
ciently constrained to test the models. Improved orientation 
constraints, perhaps from additional radio and X-ray obser- 
vations, are needed. A parallax measurement would also be 
very valuable for strengthening the use of this pulsar to test 
luminosity models. 

4.3. Other Nearby, Energetic Radio Pulsars 

The DM-distance cut selects three additional young ener- 
getic radio pulsars, but these, plotted on Figure lb, may well 
have typical LAT pulsar luminosities. First, PSR J0834— 4159 
is undetected with a flux limit giving ~ 1/3 Lj^heu at its DM 
distance. Like J0745— 5353, this pulsar's DM distance es- 
timate is dr amatically smalle r in th e CL02 model (1.6 kpc) 
than in the ITavlor & CordesI (fT99l TC93 model (9.7 kpc), 
due to inclusion of nearby HII complexes. However, un- 
like J0745— 5353, even a modest factor of 2 increase in the 
true distance would make the present flux bound unconstrain- 
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ing. Geometrical data on this pulsar are limited. It has pulse 
components separated by » 165°, so it is a likely interpulsar 
dWeltevrede and Johnston ^ 2008). However, unlike other pul- 
sars discussed in this paper it has very little linear polarization, 
so the geometry has not been confirmed by RVM modeling. 
In sum, this is not a strong case for a sub-luminous pulsar. In- 
deed, the limited radio information suggests a near-orthogonal 
rotator and thus 7 emission beamed toward Earth. Continued 
LAT exposure, and improved radio data can help clarify the 
situation. 

Statistically significant LAT flux is found in the direction 
of PSR J0857-4424. Like the other pulsars in this region, the 
DM-estimated distance had a major adjustment and should 
be considered uncertain. In addition, background systematics 
can perturb the LAT flux estimate. We conclude that this pul- 
sar is likely not sub-luminous, although a LAT pulsed detec- 
tion is required for firm conclusions. Unfortunately the lack of 
significant linear polarization in this pulsar will make it very 
difficult to extract detailed radio geometry constraints. 

The last nearby energetic object is PSR J1722— 3712 
(B1719-37) at 1.9kpc (CL02). This interpuls ar has been 
the subject of a very careful polarization study by lKeith et alJ 
(HOIQ) who find (EWOl-coiTected) angles a = 89.3 ± 0.1° 
and C, = 83.9 ± 0.3°. Thus we have high confidence that 
this is an orthogonal rotator and we expect to see visible OG 
emission. The LAT does, in fact, provide an unpulsed detec- 
tion. There are no particular issues with DM in this direction 
and the LAT detection is consistent with a point source, lo- 
calized to the radio pulsar position and having a pulsar-like 
7-ray spectrum. The inferred luminosity is quite consistent 
with Lyjieu- Thus this object has a well constrained geometry 
indicating that the outer magnetosphere 7-ray beams should 
be visible and we do indeed detect the source. Interes tingly, 
of the five d < 3 kpc interpulsars in the study of Kei th et alJ 
dloTol) three are LAT detected (PSR J1057-5226=B1055-52 
at 0.7kpc, PSR J1722-3712=B1719-37 at 1.9 kpc and PSR 
J0908 -49 1 3=B0906-49 at 2.5 kpc) while the other two (PSR 
J1549-4848 at 2.7kpc and PSR J1739-2903=B1736-29 at 
2.5 kpc) have at present high upper limits or weak detections, 
quite consistent with L^,heu- 

LAT pulse searches for these objects and similar interpul- 
sars with well constrained geometries can provide important 
model tests. Note that with ( « 90° we are probing emis- 
sion very near the null charge surface for these objects. As it 
happens, the presence of 7-ray OG emission in this region is 
sensitive to the magnetospheric currents (Hirotani 2006). The 
shape and phase of the LAT pulsations allow us to trace the 
emission to particular magnetosphere zones. When kinematic 
distances are also available, we then have the actual luminos- 
ity of these zones, a particularly powerful model constraint. 
Unfortunately such constraints are very difficult to obtain for 
the nearly-aligned pulsars which are plausibly associated with 
the sub-luminous pulsar class discussed here. 

5. CONCLUSIONS AND FUTURE PROSPECTS 

We have examined the set of 12 young, energetic radio pul- 
sars with d < 2 kpc. Most follow the standard LAT pattern 
of powerful, efficient high altitude 7-ray emission with lu- 
minosity - L^Meu- However PSR J0659H-1414, at 1/20*^^ 
of this luminosity, is a clear outlier. Present flux bounds 
make a good case that PSR J0538H-2817 is also lOx sub- 
luminous. Although they lack parallaxes, PSRs J0745— 5353 
and J1740H-1000 are also Ukely to produce < O.lL-y^heu- PSR 
J0358H-5413 has a parallax, but a less restrictive bound, at 



present. Thus we find that 4 and perhaps 5 of our nearby ener- 
getic sample may be members of this sub-luminous class, al- 
though for several of these pulsars the case could be strength- 
ened with parallax distance measurements. 

We have attempted to determine whether this may be at- 
tributed to beaming away from Earth, as would be expected 
for high-altitude (OG) emission and aligned or anti-aligned 
spin geometries. Unfortunately, precise geometrical con- 
straints are very difficult to obtain for aligned rotators. We 
do find that for our sub-luminous pulsar candidates the present 
radio constraints prefer aligned geometries. The best cases are 
probably PSR J0538H-2817 and PSR J0745-5353. Although 
it lies below the E cut for our sample set, PSR J1932H-1059 
also seems sub-luminous and aligned. However, no one case 
for alignment is compelling. 

The converse argument is in somewhat better shape. When 
we know that the pulsar is orthogonal, we seem to see 7-ray 
emission at the expected L-^^heu- PSR J1722— 3712 is an ex- 
cellent example. Although the undetected PSR J0834-4159 
may be orthogonal, its distance estimate is particularly uncer- 
tain, and its present flux limit is quite likely compatible with 

In the population synthesis of 'W atters & Romanil ( 1201 Ih it 
was concluded that for OG geometries « 30% of the radio- 
selected Log(ij) > 33.5 pulsars should be undetected in the 
7-rays, simply due to beaming. Our present sub-luminous 
pulsar fraction, (4 to 5)/12 pulsars, is consistent with this ratio, 
given the small number statistics. However we would not ex- 
pect many more sub-luminous pulsars unless the radio beams 
are substantially larger than and/or 7-ray beams are substan- 
tially smaller than assumed in these beaming computations. 
The fact that we see evidence for high altitude (/ilc > 0.1) ra- 
dio emission in several of our pulsars supports the presence of 
wide radio beams and a somewhat larger sub-luminous pulsar 
fraction; Ravi, Manchester & Hobbsl (120 10) have also argued 
that the statistics of radio and 7-ray detections imply wide ra- 
dio beams for young pulsars. 

The detection of PSR J0659H-1414 at L-y,,ie„/20 presents 
an important addition to this beaming picture. This low flux 
level, together with the 7-ray pulse's unusual phase and soft 
spectrum suggest atypical magnetospheric emission. If this 
is an aligned (or anti-aligned) rotator, this must be low alti- 
tude r < tnc emission. A set of J0659H-1414 type pulsars is 
certainly needed to probe the origin of this low flux. 

If we assume that several of our sub-luminous candidates 
join the J0659H-1414 class, we already have hints to the com- 
mon features. Certainly non-orthogonal geometries seem pre- 
ferred, although better geometrical constraints are needed for 
most sources to establish them as aligned rotators. Other plau- 
sible peculiarities for PSR J0659H-1414 exist; for example it 
has the lowest light cylinder field Slc of any LAT radio pul- 
sar Table 1 lists Blc for our candidates. No strong trend is 
evident, and certainly the few kG fields of PSRs J0358H-5413, 
J0538H-2817 and J1740H-1000 are quite typical of those of de- 
tected LAT pulsars. For these objects, at least, orientation 
seems more promising. 

Since several of our sub-luminous candidates have limits on 
luminosity approaching the J0659H-1414 level, it will be im- 
portant to see if increased LAT exposure or pulsed searches 
can detect evidence of similar low level 7-ray emission. Any 
such pulse detections for these sources must then represent 
an atypical 'sub-luminous' mechanism. The phasing of such 
pulses can be used to cement the geometrical location and 
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their spectrum and dependence on spin properties should help 
lock down the emission physics. The prospect of using the 
LAT to probe a second domain of pulsar particle acceleration, 
in addition to the established high luminosity outer magneto- 
sphere emittors, is very exciting. 
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